Quantum coherence and control in one- and two-photon optical systems 
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We investigate coherence in one- and two-photon optical systems, both theoretically and ex- 
perimentally. In the first case, we develop the density operator representing a single photon state 
subjected to a non-dissipative coupling between observed (polarization) and unobserved (frequency) 
degrees of freedom. We show that an implementation of "bang-bang" quantum control protects pho- 
ton polarization information from certain types of decoherence. In the second case, we investigate 
the existence of a "decoherence-free" subspace of the Hilbert space of two-photon polarization states 
under the action of a similar coupling. The density operator representation is developed analytically 
and solutions are obtained numerically. 



[Note: This manuscript is taken from the author's un- 
dergraduate thesis (A.B. Dartmouth College, June 2000, 
advised by Dr. Walter E. Lawrence), an experimental 
and theoretical investigation under the supervision of Dr. 
Paul G. Kwiat. 1 ! 



I. INTRODUCTION 

Decoherence in two-state quantum systems is a sig- 
nificant obstacle to the realization of proposed quan- 
tum information technologies. Coupling between quan- 
tum bit ( "qubit" ) states and unobserved environmental 
degrees of freedom leads to decoherence effects which 
limit the practical implementation of proposed quantum 
algorithms Q. Photon modes, including polarization 
and spatial modes, provide an easily accessible system 
in which simple quantum circuits can be investigated 
P, 0, m. Here, we examine the process of decoherence 
by subjecting single photons to a controllable birefrin- 
gent "environment" and observing the evolution of the 
polarization state. 

In section |J we investigate the evolution of a sin- 
gle photon state under the action of a unitary coupling 
between polarization and frequency modes. Such non- 
dissipative "phase errors" give rise to decoherence effects, 
whereby the photon evolves from a definite polarization 
state to an unpolarized state. We then introduce and 
examine an optical implementation of so-called "bang- 
bang" quantum control of decoherence by rapidly ex- 
changing the eigenstates of the coupling operation ||. 
The term "quantum control" is justified since such an op- 
eration will be shown to preserve a coherent polarization 
state in some special cases, and to reduce decoherence in 
more general cases. 



In section [II we will introduce a particular two- 



photon polarization-entangled state that, due to its sym- 
metry properties, is immune to collective decoherence 
of the type mentioned above. That is, this state is a 
decoherence-free subspace (DFS) of the Hilbert space of 
photon polarization || Photon pairs entangled in 
both polarization and frequency degrees of freedom, such 
as hyper-entangled photons produced in down-conversion 
sources (see further complicate this particular de- 

coherence mechanism . In particular, energy conserva- 
tion imposes frequency correlations which affect the co- 
herence properties of these two-photon states. In the 
experimental case, it will be shown that this frequency 
correlation can be effectively suppressed and the DFS re- 
covered by a simple and physically intuitive modification 
of the apparatus. 



II. COHERENCE AND "BANG-BANG" 
CONTROL: THE ONE-PHOTON CASE 

A. Review of single-photon decoherence 

In this section, we will show that non-dissipative (uni- 
tary) coupling between photon frequency and polariza- 
tion in a birefringent "environment" followed by a trace 
over frequency leads to decoherence. 2 These results are 
relevant to the study of decoherence in quantum systems 
which arises when coupling to environmental states, fol- 
lowed by a trace over those degrees of freedom, leads to a 
loss of phase information between qubit basis states. In 
the following argument, the term "environmental" will 
be used to describe the coupling of information-carrying 
states (photon polarization modes) to degrees of free- 
dom which are not utilized for information representa- 
tion (photon frequency modes). It will be shown that 
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2 This effect has a well-known counterpart in classical optics 
whereby quasi-monochromatic light composed of uncorrelated 
frequency components loses the ability to interfere with itself 
when polarization modes are separated beyond the coherence 
length of the incident light, as in an unbalanced Michelson inter- 
ferometer (see jij, §7.5.8). 
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even non-dissipative coupling between these modes leads 
to a loss of information in the qubit states. 

A single photon characterized by its frequency spec- 
trum and polarization can be represented by the state 
kct 



2 r 



(1) 



where \xj), j G {1: 2}, are orthonormal polarization basis 
states with complex amplitudes Cj , and A(w) is the com- 
plex amplitude corresponding to the frequency to, nor- 
malized so that 



J dw|A(w)| 2 = 1. 



(2) 



For simplicity, we assume that the frequency spectrum is 
independent of polarization, so that we need not index 
A(w) by polarization mode. In physical terms, this moans 
that polarization and frequency are not entangled. 

Since represents a pure state, the density operator 
can be written as p = (see O], §3.4), so that we 

have the initial state 



Puj{x = 0) = ^ °i c j\Xi){Xj 



dwida>2A(cJi)A*(w2)|wi)(u;2|. (3) 



The subscript indicates that p u includes frequency de- 
grees of freedom. The usual 2x2 density matrix repre- 
senting the polarization state of the photon is given by 

Pij — C{Cj. 

We will now seek the dependence of the density op- 
erator p(x) on some spatially extended "environment," 
by introducing the operator U(x) which we require to be 
both linear and unitary. Furthermore, we demand that 
U(x) have eigenkets \xj) with (frequency-dependent) 
eigenvalues Uj(u>,x) where, as before, j £ {1,2} indexes 
the polarization mode. Since U(x) exhibits frequency- 
and polarization-dependent eigenvalues, it introduces a 
coupling between qubit states (polarization modes) and 
non-qubit states (frequency modes). As shown in Ap- 



pendix Al, U(x) introduces pure "phase errors" on po- 
larization modes \Xj)i so that 



partial trace over frequency degrees of freedom: 
p{x) = { dLu(uj\XJ(x)p(0)U f (x)\Lj) 



(5a) 



= °i c *Mi){xj\ 

x / <L,|A(w)|Vfc < < w ' x >-^< & '' !e H. (5b) 



Since this process realizes the entangling of qubit states 
with environmental states (frequency modes which are 
not involved in information representation or manipula- 
tion) , we expect to observe decoherence effects in the off- 
diagonal elements of the qubit state (polarization) den- 
sity matrix ||. In order to investigate the coherence 
properties of a photon under the action of such an oper- 
ator, we observe that in a completely mixed (or "deco- 
hered" ) state, the density matrix is simply a multiple of 
the identity matrix. The diagonal elements of p(x) are 
indexed by i = j. For these values of i and j, the argu- 
ment of the exponential in Eq. ^ vanishes and by the 
normalization condition (Eq. ^), we have pu(0) = pu(x). 
In other words, the diagonal elements of p(x) are unaf- 
fected by U(ir). 

The off-diagonal elements of the density matrix are 
indexed by the values i = l,j = 2 and i = 2,j = 1. To 
compute these explicitly, we must define the functions 
<Pj(aj, x). In the case of a single photon in a birefringent, 
linearly dispersive crystal of thickness x (e.g., a quartz 
crystal) , 



<Pj\UJ, x) = L 



(6) 



where rij is the index of refraction corresponding to po- 
larization j. 3 Making the substitution r = ( -" 2 ~" 1 ^ and 
writing f((j) = |A(w)| 2 , the off-diagonal elements of the 
density matrix are given by 



(7a) 



Pu{x) = Pi2(0) / du>f{io)e' 



P2i(x) = P2i(0) / dw/(w)e l 

-' CO 

= p 2 i(0)^(r) 



(7b) 



where T{t) is the Fourier transform of f(u>) (up to a 
constant, depending on convention). 

Evidently, p^(0) = p(Q) implies that p'(x) = p(x), so 
that Eqs. [Fa] and |7^ preserve the required hermiticity of 



U0r)| Xj 



Xj) 



(4) 



where the phase factor ipj(u),x) is a real- valued function 
of uj and x. 

Now we can calculate the spatial dependence of the 
density operator under the influence of U(x) and follow- 
ing a frequency- insensitive measurement, effected by a 



Here, we neglect dispersion, the variation of n with u>. By as- 
suming \xj) to be an eigenstate of U (Eq. W), we require the 
optic axis of the birefringent element to be aligned with one of 
the polarizations \xj)- In the case where the represent cir- 
cular (elliptical) polarizations, this means that the environment 
is optically active as opposed to (as well as) birefringent. 
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p. In terms of An = ni — n\ and T(r), the density matrix 
at some later position x becomes a simple modification 
of the density matrix at x = 0: 



p(x) 



Pn(0) fta (0)^*(s|2) 



P2l(0)^(^) 



p 22 (0) 



(8) 



Eq. U allows us to calculate the density matrix repre- 
senting the polarization of an optical qubit subjected to 
a non-dissipative frequency-polarization coupling, given 
the initial state density matrix, p(0), (written in the basis 
of eigenmodes of the coupling) and the functional form 
of the frequency-amplitude function A(w). 4 In general, 
|A(o;)| 2 is peaked at some central value oj with a finite 
width Slo. Under these conditions, the Fourier transform 
T{t) (and the off-diagonal elements of p) will fall signifi- 
cantly on the time scale r c ~ In other words, there is 
no coherent phase relationship between polarization ba- 
sis states |xi) and |xa) after the photon has traveled a 
distance x > through the environment XJ(x). 

To summarize, we have developed a prescription for 
calculating the polarization state density matrix repre- 
senting a photon in a birefringent environment. From 
this expression, we see that a photon with definite polar- 
ization phase information before entering such an envi- 
ronment will, in general, lose phase information between 
polarization basis states after traveling a suitably long 
distance. The physical mechanism which destroys such 
phase information is the entanglement between polariza- 
tion modes and frequency modes, which is an instance 
of the entanglement between qubit states and environ- 
mental degrees of freedom which are not useful for infor- 
mation manipulation. Because we trace over frequency 
(environmental) degrees of freedom in this case, the phase 
coherence between the qubit basis states is destroyed. 



B. "Bang-bang" quantum control of decoherence 

In the previous section, we showed that a qubit con- 
sisting of photon polarization modes \xj) may lose its ca- 
pacity for information representation under the action of 
a unitary operator, V(x), which entangles frequency and 
polarization. In other words, decoherence occurs even 
for non-dissipative "phase errors" in which the qubit ba- 
sis states are eigenstates of the environmental coupling, 
and there is no chance for a bit flip error Q|. 5 

In ||, the authors describe this result in the general 
case of qubit states coupled to environmental degrees of 
freedom and introduce a scheme for reducing such de- 
coherence effects via rapid, periodic exchanging of the 



4 Actually, we need only know the complex square of the 
frequency-amplitude function, |A(w)| 2 . 

5 Note however, that phase errors in the basis {|xj)} mav appear 
as bit-flip errors in another basis, e.g. { ^1 X2 ^ }. 



eigenstates of the environmental coupling. In our case, 
this corresponds to an interchange of the polarization ba- 
sis states \xi) and \x2), which can be accomplished by 
an appropriate reflection or rotation operation. These 
exchanges are rapid in the sense that the period of flip- 
ping must be short compared to the decoherence time 
Tc ~ 3w an( ^ tire time scales of any other relevant de- 
coherence mechanisms (e.g. scattering, dissipation). So 
called "bang-bang" control reduces the degree of deco- 
herence by averaging out the time-dependent coupling 
between qubit states and environmental states. 

In this section, the previous results are expanded to 
include an implementation of "bang-bang" quantum con- 
trol by periodic exchanging of the environmental eigen- 
states \xj}- The mathematical formulation is consider- 
ably simplified if we move to the discrete regime in which 
the "environment" acts in a stepwise fashion, and the ex- 
change of eigenstates also occurs in discrete steps. 6 To 
this end, we define R to be a reflection (or rotation) op- 
erator which exchanges the eigenstates of the frequency- 
polarization coupling: 



R|x0 
R|X2) 



±lxi) 



(9a) 
(9b) 



Here, the plus (minus) sign refers to a reflection (rota- 
tion). The choice of sign will not affect the results of a 
measurement in the experimental case, and for concrete- 
ness, we take the sign to be positive. 

We defi ne U , a "step-wise" operator analogous to U(x) 
in Sect. II A and also define C to be the thickness of 



the birefringent element (in our experiments, usually a 
quartz crystal); as before, rij is the refractive index corre- 
sponding to polarization state j, and An = ni — n\. This 
gives a relative phase shift of ip(uj) = f2{^) — ¥?i(k>) = 
at, so that we have 7 



' CAn\ 



U\ X i) = \xi) 



U\xi) 



|X2). 



(10a) 
(10b) 



Eqs. |9a-9b and 10a - 10b have the important conse- 
quence that 



(KU) 2 = e Muj} 



_ -itp(u>) 



(WtRt 



(11a) 
(lib) 



where it is understood that such operator identities are 
only meaningful wh en the operators are applied to kets or 
bras (see Appendix A 2 a for a proof of these identities) . 



Developing the general continuous case in analogy with Sect. 



II A 



by including a rotation operator acting simultaneously with 
U(x) is considerably complicated, since this rotation operator 
will not, in general, commute with U(x) for different x values. 



Here, we omit the phase shift of ipi(uj) 
common to Eqs 



] ui, which is 



10a and lObl. Such a global phase shift is 
never observable and cannot affect the coherence properties of 
the polarization states. 
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In this discrete regime, we consider "rapid" exchang- 
ing of eigenstates to correspond to the alternating action 
of the operators U and R. Regarding U and R as optical 
elements (for example, a quartz crystal and a half- wave 
plate) , we consider a cavity scheme in which a single pho- 
ton passes through this two-element unit N times. 

Symbolically after N passes through the system, we 
have (compare Eq. |5a|), 



p(N) = J dcj(u/| (RU) N p u (0) (WtRt; 



N 



(12) 



Application of identities 11a and |llb| leads to a simpli- 
fication, depending on whether N is odd (N = 2m + 1) 
or even (N = 2m). 



p(2m) 



, 2m 



du(u\ {RU) 2m PuJ {Q) (wtRt)"" \lu) 
p(N = 0) 



(13a) 



p(2m + 1) = duj(uj\{RU) 2m+1 Pu (p)(U*Ri 



,2m+l 



dw(w|e am< *'^HWpa,(0)W t R t e- amiv ^|w> 
dw{oj\BLi Pui (0)U^\oj) 



p(N = l). 



(13b) 



Note that p(N = 1) is related to the continuous case of 
Sect. llIA| by 



P (N = 1) = /" dw(w|Rp(x = £)Rt 



(14) 



These relations lead to a simple analog of Eq. bra (where 
we use PuJ (N — 0) «-> p u (a; = 0), given by Eqf"g, as the 
input state in both cases): 

2 

p(N = 0) = J2 c i$\Xi){xA (15a) 

2 



The subscript qc (for "Quantum Control") indicates 
the case in which R is included. Note that p QC is depen- 
dent only on the parity of N , so that under the influence 
of the quantum control procedure, any number of cycles 
through this particular birefringent environment causes 
at most a dec oherence equal to that of the first pass. In 
fact, Eq. 16a shows that an even number of cycles causes 



no decoherence whatsoever. 

By making the identification x — NC (the distance 
traveled through the crystal after N pass es) in Eq. ||, we 
can dire ctly compare the results of Sect. |I A| with those 
of Sect. flBl: 



p(N) 



Pn(0) Pi 2 (0)^* (JVr ) 
p* 12 (0)F (Nt ) p 22 (0) 



(17) 



Letting r c be the coherence time such that T(t) is 
significantly reduced for r > r c , then p{N) represents 
a partially mixed state for N > On the other 

hand, Pqc{N) never loses all phase information in the 
off-diagonal elements since N — 2m recovers the initial 
state for all m. In this sense, R acts as a quantum control 
element, suppressing decoherence and maintaining state 
purity indefinitely. 

We may ask what happens if the strength of the envi- 
ronmental coupling varies with x, so that U varies on a 
time scale comparable to the period of the exchange op- 
erator, R. For simplicity, we model this "slow-flipping" 
case by introducing two distinct operators which differ 
only in their eigenvalues (not their eigenstates). In more 
concrete terms, we consider two coupling operators, U\ 
and U2, characterized by phase errors <pi{u>) and tpi{w) 
so that, for j = 1, 2: 



Uj\xi) = \xi) 

Uj\X2) = 



X2) 



(18a) 
(18b) 



The exchange operator R acts in alternation between 
each of these two operators, so that the strength of the 
coupling changes faster than the eigenstates of the cou- 
pling are exchanged. After N passes through the system, 
we have 

p(N) = f ^{w^UiU^ p u {U){U\ut) N \u) (19a) 



p(N=l) = C *S* f du\A(u)\ 2 (RU) \xi)(Xj\ (Z^R 1 ) ■ p QC {N) = /dc<c|(Ri/ 2 Ri/ 1 )^(0)(W 1 t RtW 2 t Rt) 



N 



(15b) 

Finally, writing out p{N) in matrix form, substituting 
T for the Fourier transform of f(u>) = \A(ui)\ 2 , and de- 
noting by To = the characteristic time parameter, 
we have 



(19b) 

Intuitively, we expect to see a partial reduction in the 
degree of decoherence when the quantum control proce- 
dure is implemented by introducing the exchange opera- 
tor. 



o (2m) ( Pll{0) Pl2(0) ' 

pu(0) Pt 2 (0)^(r ; 

. Pl2 (o)r* (r ) p 22 (o) 



p QC (2m + 1) 



(16a) 



(16b) 



8 Here, "mixed" refers to the fact that the off-diagonal elements 
of p are small so that Tr(p 2 ) < f . If the initial state, p(0), has 
diagonal elements of 1/2, then p represents a completely random 
ensemble for N 2a . 
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To see these results quantitatively, we make the follow- 
ing substitution which reduces this problem (mathemat- 
ically) to the previously solved one (see Appendix A 2 b 
for a derivation of these relations): 



u 



so that, 



U\ X i) 

u\xi) 



IXi> 

e t(pi(uO+v> s («0)|^ 2 \ 



U QC \xi) = 1x0 

U QC \X2) = e '(^M-^(-))| X2 ). 



(20a) 
(20b) 



(21a) 
(21b) 

(22a) 
(22b) 



With these relations, each operator is mathematically 
analogous to the operator U defined by Eq. EL wi th <p(uj) 



redefined as appropriate given relations 21af 22b and the 



details of the particular environment. In the case of birc- 
fringent crystals, the elements corresponding to U\ and 
U2 differ only in their thickness (not in their orientation 
with respect to the polarization basis states \xj))- Defin- 
ing Cj to be the thickness of a crystal represented by Uj , 
we have 



,'CjAn 



(23) 



We see that before introducing the exchange operator, 
the birefringent environment acts effectively as a single 
birefringent crystal of thickness L\ + £2 (Eq. 21b| ). Af- 
ter implementing "bang-bang" control by introducing the 
exchange operator R, the environment looks e ffect ively 
like a single crystal of thickness C\ — £2 (Eq- |22b| ). In 
the absence of quantum control, the photon is subjected 
to the sum of the individual phase errors; in the presence 
of a quantum control element, the photon is subjected 
to the difference of the individual phase errors. In this 
case, where the strength of the environmental coupling 
changes on roughly the same time scale as the exchang- 
ing of eigenstates, the quantum control procedure reduces 
but does not eliminate dccohcrcncc effects. 



Experimental Demonstration of "Bang-Bang" 
Control 



It will be convenient in our experimental analysis to 
measure the visibility, traditionally defined as 



V 



N„ 



N„ 



N„ 



N„ 



(24) 



where N max (N m i n ) denotes the maximum (minimum) 
number of photons passing through orthogonal settings 
of a polarizer. Denoting these polarization analyzer set- 
tings by I A) and |A), we define visibility at the single 



PBS 



Recycling Mirror 




Compensator 



C\ " *" Polarizer 



Detector 



FIG. 1: Schematic diagram of experimental apparatus for pro- 
ducing controllable decoherence in the H/V basis. PBS indi- 
cates a Polarizing Beam Splitter; R indicates the (removable) 
exchange element, a quarter- wave plate in this case. The in- 
put state is |45°) (= ]H) ^ V) ). 



quantum level by calculating probabilities from the den- 
sity operator, p(x) 9 : 



V A (a;) = |<A|p(x)|A> - <A|p(a;)|A>|. 



(25) 



Note that we no longer include the denominator corre- 
sponding to the total number of photons in Eq. E4, since 
the denominator here is the trace of p in the {[A), |A)} 
basis, which is unity. V\(x) quantifies the similarity be- 
tween the state represented by p{x) and the state |A) 
in the following sense: if p(x) represents an unpolar- 
ized state, it is a multiple of the identity matrix, and 
V\(x) = 0; at the opposite extreme, if the photon is 
in a pure state with polarization A, p(x) — |A)(A|, and 
V\{x) = 1. If we choose to measure with respect to the 



state I A) =^{\Xi) 
for a proof) 



IX2)), we find that (See Appendix 



V(x) = 2Re((xMx)\ X 2)) 



(26) 



where we have dropped the subscript A for this choice 
of measurement, and Re indicates the real part. In 
this form, it is clear that visibility V(x) depends on the 
off-diagonal elements of p{x), and so corresponds with 
our previous notion of coherence. If we write p in the 
H/V basis, measuring the visibility with respect to |A) = 
^(\H) + \V)) = |45°> and |A) = ±(\H)-\V)) = |-45°) 
gives us information about the off-diagonal elements of 
P- 

To see the results of the previous sections experimen- 
tally, we seek a systematic and controllable implementa- 
tion o f the operators U and R in the cavity scheme of 
Sect. II B. For the first case, in which there is a single 



operation which produces a phase error, we use a Michel- 
son polarization interferometer as a tunable decoherence 



9 The notation {|A), |A)}, is intended to suggest that these polar- 
ization modes may, in general, represent any elliptical orthonor- 
mal basis states. 
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mechanism in the H/V basis. A polarizing beam-splitter 
(PBS) separates the path of V- and ii-polarized light. 
The unbalanced arms of the interferometer produce a 
controllable, ^-dependent phase shift between \H) and 
\V) (see Fig. |). 

A photon in cither path is reflected back to the PBS 
and retraces its path to the input. The difference in path 
lengths between arms of the interferometer will be de- 
noted by x. The unbalanced interfero mete r is the optical 
real ization of the operator U of Sect. 
II B. 



II A 



and U of Sect. 



An R = 0.97 reflector at the input recycles the light so 
that a photon passes through the unbalanced arms of the 
interferometer more than once, in general. 10 A partially 
reflective out-coupling mirror directs a photon out of the 
interferometer with a probability of R — 0.04 at each 
pass. 11 

We send in light from a greatly attenuated diode laser 
at A = 670 nm, pulsed at 100 kHz. In our experiment, 
we observe photons passing through the interferometer 
< 10 times. The total distance traveled by a photon 
passing through our interferometer 10 times is ~ 5 m, so 
that any photon spends at most ~ 10~ 8 s in the system. 
EG&G SPCM-AQ silicon avalanche photodiodes (dark 
count < 400 s _1 , efficiency ~ 60%) are used as photon 
counters. Using electronic timing techniques, we record 
photon counts only in a coincidence window of 50 ns be- 
tween the pulse generator and detection event. These 
timing techniques, together with spectral filters at the 
detectors, reduce background coincidence counts to < 1 
s _1 . We graphically display photon counts vs. arrival 
time so that photons leaving the system after N cycles 
are recorded in the Nth peak. In a typical experiment, 
our photon counters register no more than 5000 coinci- 
dence counts per second, and we conclude that ~ 10, 000 
photons per second leave the system (~ 0.1 per pulse). 12 
At the output, we count for a few seconds (less than 
a minute in most cases). Although we cannot determine 
the unknown polarization state of a single photon, we can 
accurately find the statistical distribution of the ensemble 
over the observation time, from which useful information 
about the density matrix can be inferred. 

By placing a A/4 waveplate with optic axis at 45° at 
the recycling reflector, we realize an optical "flipping" 
operator. The A/4-reflector-A/4 combination takes \H) 
to \V) and vice versa. This is the optical realization of 



the exchange operator, R in Sect. II B. 

By recording photon counts in coincidence with the 
diode pulse generator and using a polarizer at the detec- 
tor, we can measure the visibility at |45°) = -^(\H) + 

\V)). The frequency spectrum of laser light is restricted 
using an interference filter with bandwidth SX = 10 nm 
and central wavelength A Q = 670 nm. 13 For our calcu- 
lations, we use a Gaussian amplitude function with fre- 
quency spread 5u> corresponding to 5X — 10 nm, central 
frequency loq corresponding to A Q = 670 nm, and nor- 
malization factor A a : 1A 



A(lj) = A a exp 



Sw/y/2 



(27) 



The corresponding coherence length of the photons is 

A 2 

given by L c ~ -gf = 45/xm, so we expect substantial 
decoherence for a path difference of ~ 4.5 /xm at cycle 
N = 10. 

W ith this setup, we can measure V vs. N as in Sect. 
|II B| . Experimental and theoretical results are plotted in 
Fig. H for the case where the arms are unbalanced by 
x = 0.00, 2.49, and 4.99^m. The experimental results 
show clearly that the introduction of the quantum control 
element results in a significant reduction of decoherence 
in good agreement with the theory. However, practical 
difficulties with the interferometer limit the quantitative 
agreement between experimental and theoretical values 
of V QC , especially at large N. We see in Fig. |2|(a) that 
the relative path difference is not after N = 10 passes 
through the interferometer due to imperfect balancing. 
By determining the coherence length from the data in 
Fig. ||, we infer that the actual path length difference in 
the interferometer was ~ ±0.7 /jm (about one wavelength 
of incident light) displaced from the nominal x- values. 
This discrepancy then determines the bounding value of 
V(N = 10) for all other plots. Theoretically, of course, 
this value is exactly unity. 

In order to model the "slow-flipping" case, we intro- 
duce an additional operator, U2, which will induce a sec- 
ond phase shift, also in the H/V basis. To accomplish 
this task, we include a birefringent quartz crystal whose 
optic axis is aligned with the axes of the PBS of the 
original setup. We denote the optical path difference in 
this crystal by xq. 15 We now use a A/2-waveplate as 



10 We use R to denote the probability of reflection at a surface. 

11 For normal incidence, R = 0.04 at each air-glass interface. 
Here, one surface of the out-coupling mirror is broad-band anti- 
reflection (BBAR) coated so that we have effectively only a single 
air-glass interface. 

12 The probability that there are two photons in the system si- 
multaneously is given by the Poisson distribution, P(2) = 

e -0.1 ( Q j\2 

2 I ' ~ 0.005. These photon numbers are comparable to 

those used in current demonstrations of free-space quantum cryp- 
tographic key distributions |l2|. 



3 Since the diode laser is weakly driven just above its lasing thresh- 
old, the intrinsic frequency spectrum of the photons is wide 
enough that these filters determine the overall frequency distri- 
bution. 

4 In practice, we find that the frequency spectrum is not as smooth 
as a Gaussian. From the visibility curves, we infer that the fre- 
quency spectrum has a substantial monochromatic (compared to 
8b>) component at u> . 

5 Instead of a single quartz crystal, we use a 1.046 mm crystal 
at 0° and a 0.850 mm crystal perpendicular to the first. The 
combination of these two crystals gives an effective path length 
such that nearly complete decoherence occurs over 10 cycles. 
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FIG. 2: Experimental and theoretical curves showing V vs. TV 
for the case where (a) x = 0.0 fj,m,(b) x — 2.49 /im, and (c) 
x = 4.99 /im, for the setup in Fig. Fy. "QC" indicates the case 
in which the quantum control procedure was implemented. In 
(a) (x — 0), the theoretical curves for both cases are constant 
at unity, since we expect no decoherence in this case. In 
(c), the observed visibility curve (without quantum control) 
falls to V ~ 0.15 on the length scale of our observations. 
We suspect that V falls to on the length scale determined 
by the monochromatic component of the frequency spectrum 
mentioned in footnote The experimentally observed values 
in the quantum control case alternate in the opposite sense to 
the theoretical values. This discrepancy is not yet understood. 



the quantum control element which exchanges the eigen- 
states, H and V, at each pass both to the left and to the 
right (see Fig. ||). For the theoretical values, we used the 
same prescriptio n as abo ve, but with effective operators 
defined by Eqs. |21a - [22b . Experimental and theoretical 
results are displayed in Fig. [|. 
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FIG. 3: Schematic diagram of experimental apparatus for the 
"slow-flipping" case. As before, R indicates the exchange ele- 
ment, a half-wave plate in this case. The birefringent crystal, 
a quartz plate, introduces a second phase error in the same 
basis as the unbalanced interferometer arms. 
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FIG. 4: Experimental and theoretical curves showing V vs. 
N for the "slow-flipping" case where xq = 1.77 /im and (a) 
xi = 2.1 fim,(b) X2 = 3.16 fim, and (c) £3 = 4.2 /im for the 
setup in Fig. 
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FIG. 6: Schematic diagram of experimental apparatus for the 
case in which the eigenstates of the polarization-frequency 
coupling can be adjusted. As before, R indicates the exchange 
element, an optically active 90° quartz rotator in this case. 
The two birefringent crystals introduce different phases errors 
whose eigenstates can be adjusted independently by changing 
their orientation. 



and 



FIG. 5: Measured (reciprocal) optical path difference inferred 
from the coherence leng ths in Fig. ^ versus expected values 
according to Sect. O H . Error bars in both directions are 
within the points. See the text for an explanation. 



In Fig. ||, we have inferred the effective rate of decoher- 
ence (a reciprocal length, in this case) in the system by 
fitting curves to each of Fig. || and (reciprocally) plotted 
these versus the rate of decoherence for a single pass pre- 
dicted by Eqs. 21t and|22bl We expect to see an effective 



decoherence rate given by the sum of xq, the decoher- 
ence rate due to the quartz crystal, and Xi, i € {1,2,3}, 
the decoherence rate due to the interferometer, when the 
quantum control procedure is not implemented; further- 
more, we expect to see a decoherence rate determined 
by the difference of xq and x when the quantum con- 
trol procedure is implemented. Here, the decoherence 
"rate" of an optical element is directly proportional to 
the reciprocal of the optical path difference between e- 
and o-polarizations. The linearity of this plot indicates 
that the q uantum control element acts as predicted in 



Sect. II B 



We are limited in the type of decoherence we can ob- 
serve by measuring V. For example, a photon in the pure 
state | H) has a visibility of with respect to the 45° ba- 
sis. We would like a more general method of quantifying 
decoherence effects. To this end, we can extend the var- 
ious methods of polarization analysis in classical optics 
to our more fundamental quantum mechanical setting. 
With this idea in mind, we define the degree of polariza- 
tion V, and the Stokes parameters (sj) in terms of p(x) 



as 



16 



^V^i) 2 + <s 2 ) 2 + (s 3 > : 



(28) 





= 2{H\p{x)\H)-l 


(29a) 




= 2(45°|p(a;)|45 ) - 1 






= (H\p(x)\V) + (V\p(x)\H) 


(29b) 


s 3 ) 


= 2(R\p(x)\R)-l 






= i((H\p(x)\V) - (V\p(x)\H)). 


(29c) 



As usual, \H) and \ V) denote horizontal and vertical po- 
larization, |45°) = 4^{\H) + \V)) denotes 45° polariza- 
tion, and \R) = ^(\H) + i\V)) denotes right-circular 
polarization. 

Note that photon counting in the appropriate polar- 
ization basis allows experimental determination of the 
Stokes parameters, ( sj). By photon counting and us- 
ing the relations 29a- [29c , the density matrix of an arbi- 
trary polarization state can be reconstructed from these 
data. 1 In contrast to measurements of V with respect to 
some particular orthonormal basis {|A), |A)}, this tomo- 
graphic technique allows complete determination of the 
polarization state. In fact, the visibility V with respect 
to any basis |A) can be found from the density matrix if 
the transformation from \xj) to |A) is known. 

In the previous cases, we used an appropriately ori- 
ented waveplate as a realization of the exchange operator, 
R. However, the choice of orientation in both cases re- 
quired prior knowledge of the eigenstates of the coupling 
between frequency and polarization. Note, however, that 
if we know only that the eigenstates of the coupling are 
orthogonal, linear polarization states, then a rotation by 
90° (as opposed to reflection about a specific axis) meets 
the requirements of an exchange operator (Eqs. ^a-^b) in 
all cases. We expect that an optically active 90° quartz 
rotator acts as a quantum control element, independent 
of the orientation of the quartz crystal that causes deco- 
herence (since the eigenstates of our birefringent quartz 



These are a generalization of the Stokes parameters of classical 
optics to the quantum mechanical case. In classical optics, these 
data can be used to reconstruct the coherency matrix, J, which 
completely characterizes the polarization state of the light (See 
@> §10-9). 



An extension of this technique of quantum tomography has been 
de mon strated by reconstruction of two-photon density matrices 
in JL3| |. This technique will be useful in the two- photon state 
discussions of Sect. III. 
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FIG. 7: Theoretical and Experimental plots of V vs. N for 
environmental eigenstates of —10° in a linear cavity (see Fig. 
^|). The input polarization is 35°. Clearly, introduction of the 
quartz rotator results in a preservation of photon polarization 
information for this choice of environmental eigenstates. 



FIG. 8: V vs. N for the case in which the strength and the 
eigenstates of the polarization-frequency coupling change. In 
the case where we include the quantum control element, a 90° 
quartz rotator, the degree of polarization is preserved. 



crystal are orthogonal and linear, i.e., along the optic 
axes of the crystal). 18 

To investigate decoherence in an "adjustable" linear 
basis, we use a cavity with a photon passing through 
1.046 mm and 0.850 mm birefringent quartz crystals at 
each cycle. Since this thickness of quartz results in sub- 
stantially faster decoherence than the unbalanced inter- 
ferometer of the previous cases, we use a narrow band- 
width filter of <5A = 1.5 nm to increase the coherence 
length of the photons. The exchange element, a 90° 
quartz rotatory is placed between the birefringent crys- 
tals (see Fig. g). At the output, we measure the Stokes 
parameters by photon counting in the appropriate bases 
and calculate V to quantify the degree of coherence. By 
orienting both birefringent crystals so that the optic axis 
makes an angle of —10° with the horizontal, and sending 
linearly polarized light at 35° into the system (for max- 
imum decoherence), we see that the rotator does indeed 
act as a quantum control element (Fig. 0). Note however, 
that we need not rotate the optically active element by a 
corresponding 10°, because its functionality is indepen- 
dent of such a rotation. Put differently, the rotator acts 
as a quantum control operator for (slow) decoherence in 
all linear bases. 

For the theoretical curves in Fig. 0, we need not modify 
the previous calculations except to take the lower (nega- 
tive) sign in the definition of R, Eq. |9b| and include the 
narrower frequency spectrum. The invariance of physical 
systems under rotation allows us to rotate the entire ex- 
perimental apparatus (the input polarization and quartz 
crystals) by 10° without adjusting the calculations. 



If the eigenstates of the environmental coupling are left- and 
right-circular polarizations, a A/2 waveplate at any orientation 
acts as the exchange operator R. 



As a further demonstration of "bang-bang" quan- 
tum control of decoherence, we seek a scheme in which 
the strength and the eigenstates of the polarization- 
frequency coupling change between exchange operations. 
To reproduce these conditions, we simply place the 
quartz crystals of Fig. |6| at slightly different angles so 
that the eigenstates of the coupling are not constant at 
each pass. The first (1.0 mm) crystal is oriented at 10° to 
the horizontal while the second (0.8 mm) crystal is placed 
at 0° degrees. To solve for V starting from the density 
operator representation, we rewrite the first crystal oper- 
ator as a matrix in the H/V basis using the appropriate 
rotation matrices. The resulting output density matri- 
ces are calculated numerically and V is extracted. The 
theoretical and experimental results are displayed in Fig. 

1- 

As a final demonstration of bang-bang control, we as- 
sert that the same technique of exchanging orthogonal 
polarization states will also compensate for decoherence 
via dissipative (i.e., non-unitary) errors. Experimentally, 
we accomplish this task by inserting a neutral-density (u- 
independent) filter into the transmission {\H)) arm of the 
balanced Michelson interferometer of Fig. [l]. This filter 
results in 35% loss (after two passes) per cycle. Denoting 
this operator as T, letting T = 1 — 0.35 = 0.65, and using 
the input state ^{\H) + 1^0), we have, after N cycles, 

a\H)+(3\V) t" | aT N \H)+(3\V). (30) 

In the case when the quantum control procedure is im- 
plemented, T and R alternate so that the state evolves 
in the following steps 

a\H)+(3\V) T , aT\H)+l3\V) 

R , aT\V)+(3\H) 

T , aT\V) + f3T\H) 

R , aT\H)+/3T\V) (31) 
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FIG. 9: Experimental and theoretical results displaying visi- 
bility vs. cycle for T = 65% transmission in the \H) arm of 
the interferometer (see Fig. [j]). The experimental data have 
large error bars since photon numbers are quite low after 10 
passes, and counting statistics (which scale as the square root 
of the photon count) become significant. As before, we expect 
that the imperfectly balanced interferometer arms will reduce 
the visibility below the theoretical values. Nevertheless, the 
principle of the quantum control for reducing the effect of 
dissipation is clearly demonstrated. 



SO that 

(TR) 2 (a|iT>+/3|V» - T(a\H) + 0\V)). (32) 

In this way, we prove the operator identity (TR) 2 = T so 
that, after an even number of cycles, the system evolves 
back to its original state but with a reduced probability of 
detection. Experimental and theoretical results for this 
case are displayed in Fig. ||. 

In summary, we have shown that coupling information- 
carrying (polarization) degrees of freedom to unobserved 
(frequency) degrees of freedom results in a loss of polar- 
ization phase information. We have shown in this sec- 
tion that an optical realization of "bang-bang" control 
of decoherence does indeed preserve photon polarization 
in the presence of such an "environment." Our imple- 
mentation of bang-bang control requires some knowledge 
of the eigenstates of the polarization-frequency coupling. 
We have demonstrated such control when the eigenstates 
of the coupling are linear and orthogonal. 19 In addi- 
tion to preserving polarization phase information in the 
presence of a frequency-polarization coupling, we have 
shown that such bang-bang control similarly reduces de- 



19 Tri f act, in the Poincare Sphere representation of polarization (See 
[JlCt] §1.4), it is sufficient to know only the plane in which the 
eigenstates of the coupling lie: a 90° rotation about the normal 
to that plane satisfies the conditions of an exchange operator R. 
The technique fails, however, for a decoherence mechanism whose 
eigenstates coincide with those of R. (i.e., which are parallel to 
this normal direction). 



coherence when the strength and/or the eigenstates of 
the coupling vary between exchanges. 



III. DECOHERENCE-FREE SUBSPACES OF 
TWO PHOTONS 

A. Polarization-frequency coupling in correlated 
photon pairs 

It has been shown that the antisymmetric 2-qubit state 
given by (in our optical notation) \ip~) = -^(|Xi)IX2) — 

IX2) |xi)) is immune to collective decoherence, i.e. deco- 
herence arising from an operation that is invariant under 
qubit permutations Q. That is, due to its symmetry 
properties, I?/ 1- ) is a decoherence-free subspace (DFS) of 
2 qubits. Therefore, it is reasonable that we expect the 
two- photon polarization state -j^(\H)\V) — \V)\H]) to be 
immune to decoherence of the type discussed in Sect. D, 
which arises from coupling to frequency degrees of free- 
dom that are unused for information representation. 

In this section, in order to investigate the possibility 
of a DFS, we will consider polarization phase informa- 
tion in a two-photon system. In non-linear optical ma- 
terials, spontaneous parametric down-conversion events 
can be used to produce spatially separated, polarization- 
entangled photon pairs with high fringe visibility J8|, |). 
Here, we extend the results of Sect, [ffl by presenting a 
method for calculating the density matrix representing 
the polarization of such a two-photon state in the pres- 
ence of a birefringent environment. Experimentally, we 
use the source detailed in || which produces correlated 
photon pairs in the two paths labeled L and R (See Fig. 

p- E>- 

In order to keep track of frequency, spatial mode (L or 
R), and polarization, we define a shorthand notation for 
indexing the 4x4 two-photon density matrix p(x). The 
element \xij)(Xkl\ is defined to be \xi) l\Xj) R{Xk\h{xi\B. 
where the subscripts L and R refer to the left and right 
photon paths respectively and |xi) and \x2) refer to the 
polarization basis states. As a concrete example, in the 
H/V basis, the state \V)l\H)r is written as \xix) an d 
\V) l\H) R {H\ R {H\ L is written as |x2i)(xn|- In this way, 
each of the 16 entries of p is indexed by four numbers, 
(ij, kl), where i, j, k, I £ {1, 2}. 

In order to satisfy energy conservation, the frequen- 
cies of photon pairs produced in a down-conversion event 
must sum to the frequency of the pump beam, and their 
resultant energies become entangled (see Jl4| §22.4). De- 
noting the pump frequency by ix> and the frequency of 
daughter photon j by loj (j G {1, 2}), we therefore require 



Mo 
~2 



uJ-2 



(33a) 
(33b) 
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so that we have u>\ + u> 2 = u in all cases. 20 This fre- 
quency "anti-correlation" will be important in an exper- 
imental demonstration of a DFS. 

Taking to be the (in general, complex) amplitude 
of the state \xij) an d denoting the frequency spectrum 
by A(w), we write the initial (pure) state at x = as 21 



where we must be careful to treat u> and i 
of the variables of integration e and e. 



as functions 



In order to represent a birefringent crystal in each path, 
we define the operator 



1-3/(0)) = ctflxy) 

de A( T + e)A(^ - e)|^ + e)|^ - e). 

(34) 

For completeness, and to make this indexing scheme 
explicit, we write the reduced 4x4 density matrix repre- 

\XijY- 



(35) 



Note that there are six independent parameters cor- 
responding to the four complex amplitudes of the 
states \Xij) constrained by the normalization condition 
\°ij\ 2 = 1 an d the freedom to choose the global (over- 
all) phase. 22 

Introducing the shorthand notation uj^ — 



senting 


the pure polarization state |\1/) = 


-T 2 






^|cn| 2 cnc* 21 cnc* 12 


c\icl 2 \ 


P = 


l*X*l = 


C2ic*n \c 2 i\ 2 c 21 c* 12 


C21C22 


Cl2C* x CxiC* 21 \c 12 \ 2 


C12C22 






\C22Cu c 22 c 21 c 22 c* 12 


\C22?) 



± e and Q 



, ^ M z (e) = ^ ± e, the density operator 

for the initial state, including frequency modes, can be 
written as 

Pu (o) = |*(o)X*(o)| 

CijC* k i\Xij){Xkl\ 



E 



x | W +)| W ->(«+ 

We have used the notation of 



de de {A(w + )A(cj-)A*( t I)+)A*(w-) 

L 



(36) 

where such frequency entan- 
glement is shown to cause nonlocal cancellation of dispersion to 
second order in optical systems. 
21 We impose the normalization condition 



U(x L ,x R ) = V L (x L )®U R (x R ) (37) 

which depends on the distance, xl [x r ), traveled through 
such a crystal in the L (R) photon path. It should be 
noted that the (£> in Eq. p9 represents a tensor product 
between the Hilbert spaces of photon polarization and 
photon frequency, whereas the ® in Eq. |37] represents 
the tensor product of photon spatial modes L and R. In 
analogy with Sect. II A, XJl(xl) associates a frequency- 
dependent phase, tfjjw, xl), to each polarization state 
j € {1, 2} in the path labeled L: 



UL(x L )\x j )L=e i ^<^\ Xj ) 



(38) 



•*)|- 



1 



A similar relation holds for the path labeled R. The func- 
tions (pj(aj,XL/ R ) characterize the phase shift induced in 
either path (L or R). In this way, we require that the cou- 
pling between frequency and polarization share the same 
dependence on the parameters x and cu in either photon 
path. This requirement is reasonable since we wish to ex- 
plore the robustness of the antisymmetric state against 
collective decoherence. 

Making a frequency-insensitive measurement of polar- 
ization after the photon in path L (R) has traveled a 
distance x^ (x R ), we have 

so that initially (*|*) = 1 and p 2 = p. 
22 For a partially mixed state, there are 15 independent parame- 
ters corresponding to 4 real diagonal elements reduced to 3 by 
the normalization constraint, and 12 complex amplitudes in the 
upper right entries. 



p{x L ,x R ) = J J du 3 du 4 (u 3 \(u i \U(x L ,x R )p LJ (0)lJ* l (x L ,x R )\uj 3 )\uj4) (39) 

which gives 

p(x L ,x R ) = J2 dAilXi^iXkil Jde |A( w +)| 2 |A( W -)| 2 e<^^ + ^)+^(-"^)-^(- + ^)-^(-"^)). (40) 

i,j,k,l 

1 



Eq. ^ gives a prescription for calculating the den- sity matrix representing polarization degrees of freedom 



12 



in the pres ence of non-dissipative "phase errors." As in 



Sect. II B , we expect to observe decoherence effects in 
the off-diagonal elements of p. Note that for the diag- 
onal elements of p, indexed by i — k,j — I (see Eq. 
|35| ), the argument of the exponential is identically zero 
so that the diagonal elements do not change. Of course, 
this result holds only for the case of a unitary frequency- 
polarization coupling in this particular orientation (i.e., 
with eigenstates \H) and \V)). 

Let us now choose the functions ipj (uj, x) to be identical 
to those of Eq. ^. In other words, we choose to realize 
a unitary polarization-frequency coupling by placing a 
birefringent crystal across both photon paths with e- and 
o-axes aligned with and \X2}- 23 Letting C denote 
the thickness of the crystal, the phase difference between 
these functions evaluated at xl = xr — £ is given by 



ip 2 (u},£) - (fii(uJ,C) = UJ 



CAn 



(41) 



In Appendix A4, the 4x4 density matrix represen- 
tation is written for the evolution of the Bell states 
|0±> = j-(\HH)±\W)) and |^±) = ±{\HV)±\VH)). 

From these, we see that the states \tp ) become "inco- 
herent" in the presence of such a frequency-polarization 
coupling. Here, we do not make a quantitative mea- 
sure of coherence (or mixture) but simply observe that 
the off-diagonal elements of the density matrix approach 
zero for C larger than the coherence length of the down- 
converted photons, given by ~ where 5oj is the width 
of the frequency spectrum |j4(w)| 2 . Therefore, these 
states lose phase coherence information in the presence 
of polarization-frequency coupling of this type: identi- 
cal birefringent crystals in the same orientation in both 
photon paths. 

On the other hand, the states l^) do not undergo 
decoherence for this particular coupling: they retain a 
definite phase relationship in the off-diagonal matrix el- 
ements. We see from the density matrix representation 
that the state -^(\HH) ± \VV)) at the input looks like 

-^(\HH) ± e^^lVV)) at the analyzer. 24 For this 
particular arrangement, a birefringent crystal with axes 
along H and V, the states \HH) ± \VV) do not deco- 
here: there is no decline in magnitude of the off-diagonal 
elements of the density matrix. 

Note, however, that neither of the states |^ ± ) can be 
considered a decoherence- free subspace (DFS). In gen- 



23 We can realize this situation experimentally either by placing 
the same crystal across both paths, or equivalently by placing 
"identical" crystals in each path and at the same orientation. 
For future convenience, we choose the second option. 

24 The phase factor arises from the pump frequency uj which we 
have assumed to be effectively monochromatic so that only a 
definite phase shift, and not decoherence, occurs. Including a 
frequency spectrum of finite width for the pump would result 
in additional decoherence effects on a different length scale, the 
coherence length of the pump beam. 



eral, with respect to another (possibly elliptical) basis 
{|A>|A», 

^=(\HH) ± \VV) * -L(|AA) ± |AA). (42) 

In particular, using the basis of left- and right-circular 
polarization, it can be shown that -^(\HH) + \ VV)) = 

±(\LR) + \RL}). So, just as the state -±(\HV) + \VH)) 
loses phase information in a birefringent crystal whose 
eigenmodes (for photons in either path) are \H) and |V), 
the state -^=(\LR) + \RL)) loses phase information in 
a crystal whose eigenmodes are \L) and \R). On these 
grounds, the only candidate for a DFS is the singlet state 
\f'~) = -j=(\HV) — \VH)), since, in analogy with the 

spin-i singlet state i(| ||) — I IT)), this state is rota- 
tionally invariant and looks the same when written in any 
basis (see jy], §3.9). However, we have already shown 
that the state \ip~) loses phase information in this sys- 
tem. 



B. Recovering a decoherence- free subspace 

The apparent contradiction between the prediction of 
a DFS and the contrary result in the previous section 
arises as a result of the frequency entanglement (anti- 
correlation) in down-converted photon pairs. We ex- 
pected to observe a DFS of two-photon states subjected 
to identical phase errors. Frequency anti-correlation 
breaks this qubit permutation symmetry since the phase 
errors are frequency-dependent. There is, however, a sim- 
ple experimental alteration which compensates for this 
frequency anti-correlation and allows us to observe a DFS 
as was originally expected. 

Global phase shifts are never observable: only the 
phase difference induced between eigenstates of the cou- 
pling operator enters into the calculations. As a result 
of frequency anti-correlation, the phase shift experienced 
by a photon in path L with frequency ^f- + e corresponds 
to that of an (identically polarized) photon of frequency 
— e in path R. 25 We have already seen that cancella- 
tion of phase terms in the exponential of Eq. ^ accounts 
for the coherence properties of a particular polarization 
state. In particular, an off-diagonal matrix element goes 
to zero if not all e terms cancel out, and the Fourier 
integral does not reduce to unity. Since frequency anti- 
correlation directly affects the sign of these phase terms, 
it also affects their cancellation or non-cancellation and 
consequently alters the predicted decoherence-free nature 
of the singlet state. 



25 This equivalence holds only in the limit that we neglect dispersion 
effects in the crystals so that An of Eq. ^ is independent of 
lu. This is true to a very good approximation over the 5 nm 
bandwidth of the photons in our experiment. 
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By inducing a second anti-correlation, which we shall 
term "path anti-correlation," we recover the proposed 
DFS. This operation consists of forcing the phase dif- 
ference in path L to be the negative of the correspond- 
ing phase difference in path R (for identical frequencies). 
By twice reversing the sign of corresponding phase shifts 
(once with frequency anti-correlation and a second time 
with path anti-correlation) we effectively produce identi- 
cal phase errors in each path, and restore the permutation 
symmetry of the polarization-frequency coupling. 

We are now faced with the experimental task of forc- 
ing phase shifts for corresponding polarizations to be of 
opposite sign in opposite photon paths. Experimentally, 
we may simply rotate the crystal in one arm by 90°. This 
operation exchanges n e and n and reverses the sign of 
An in Eq. 



41 



2(1 



To this end, we introduce a new (but familiar looking) 
operator defined by 



tj(x L ,x R ) 



XL) 



Ur(x R ) 



(43) 



\J-^(xl) represents the crystal in path L which has been 
rotated and is now characterized by the phase function 
(f>j so that 



Ufi|Xi)fl 



Xj)L 

X])r- 



(44a) 
(44b) 



The phase functions <p% and (pi in path L are related to 
the phase shifts ipi and <f2 in path R by 



(45) 



(compare Eq. ^T|). Keeping in mind that photon states 
with spatial mode L are indexed by i and k, we compute 
the analog of Eq. |4^. 

p(x L ,x R ) = 

E wUXii^Xki] I de {|A( W +)| 2 |A(, 

A 1,. ( •* 



i,j,k,l 



x e 



*»}. 



(46) 

The 4x4 matrices representing the evol ution of the 
four Bell states are displayed in Appendix A4. From 
these, we see that the states at the input become 
mixed, or incoherent, in that the magnitude of the off- 
diagonal elements approaches zero. On the other hand, 



26 In general, we can produce a similar frequency-polarization cou- 
pling with any elliptical eigenmodes A), |A) by rotating these 
eigenstates into \H) and \V) by passive optical elements (a uni- 
tary transformation), introducing a birefringent crystal aligned 
along H and V, and subsequently inverting the polarization 
transformation with more optical elements. Also, in this more 
general case, phase anti-correlation is realized by rotating the 
birefringent crystal in one arm by 90° . 



Down-Conversion Crystal (BBO) 
L- 



Birefringent Crystals 



Pump (UV) 




Phase Tuner 



Correlated Pairs 



Analyzers 



Detectors 



FIG. 10: Schematic diagram of an experimental apparatus for 
producing and analyzing correlated photon pairs. The phase 
tuner is used to produce the desired input state. The bire- 
fringent crystals have identical thickness, but can be oriented 
independently. For the experimental results, the crystal in 
path L was rotated by 0l = —17°, and the crystal in path R 
was rotated by 8r = —107°. The birefringent crystals give an 
effective optical path difference between H and V of ~ 140A o , 
where X is the central wavelength of the down-converted pho- 
tons (702 nm, in our case). 



the states |V' ± ) are possibly transformed at the pump 
frequency, but do not decohere. In particular, for the 
initial state p(0) = \if)~)(ip~\, the output state is 



P(0) 



A> 
o 












0/ 



(47) 



The state \ip~) = ^=(\HV) - \VH)) at the input looks 
like ^(|-ffV) - e l ^^\VH)) at the output. 27 



Recall now that the singlet state, \ip ) can be written 
as -^=(|AA) — | AA)) where |A) and |A) are the (in general, 
elliptical) eigenmodes of any birefringent environment. 
We conclude, therefore, that the state \4>~) is generally 
immune to decoherence of this type (i.e., collective deco- 
herence, produced here with perpendicular birefringent 
elements in each arm). By imposing an experimental re- 
quirement which restores the symmetry of the coupling 
with respect to photon path permutation, we have recov- 
ered the singlet state as a DFS with respect to collective 
decoherence. 

The technique of quantum tomography allows exper- 
imental reconstruction of two-photon polarization state 
density matrices by an appropriate choice of polarization 
measurements (on members of an ensemble) |T^| . With 
the experimental setup shown in Fig. [l^ and using this 
technique, we compare the theoretical predictions with 
experimentally observed density matrices. Details of the 
experiment and results will be reported elsewhere [ JL6| . 
Here, we simply present the theroetical predictions. 



27 In our experiment, we adjusted C (by slightly tilting one of the 
quartz elements) such that this phase factor is (modulo 2tt). 
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Input Output 




FIG. 11: Real part of polarization state density matrices for 
Bell states at the input of the experimental setup of Fig. [jlj 
as predicted by Eq. Imaginary parts are zero for all 

elements. 

We seek to quantify decoherence properties in a single 
parameter. To this end, we use the Fidelity of the trans- 
mission process, defined here as F — Tr(p out pi n ) where 
Pin (pout) represents the polarization density matrix at 
the input (output) of the system fl^| . For a mixed input 
state, we must use F = [Tr(y/^/p~p outy /p~\ Q. A 
Fidelity of 1 indicates a decoherence-free process. Nu- 
merical evaluation of Eq. |^| allows us to calculate the 
Fidelity of the transmission process for coupling eigen- 
states rotated by an arbitrary angle. Such an analysis 
reveals a Fidelity of 1 for |?/> - ) in a birefringent environ- 
ment with arbitray eigenstates, confirming the prediction 
that the singlet state is a decoherence-free subspace. 

In Fig. O, we present a graphical representation of 
the density matrices for each of the four Bell States pass- 
ing through birefrin gen t crystals in each photon path as 
calculated from Eq. Eq. 

We see that the singlet state does indeed maintain 
full coherence in the presence of a unitary frequency- 



polarization coupling. Since polarization degrees of free- 
dom are analyzed and frequency degrees of freedom are 
unobserved, such a result demonstrates the robustness of 
a decoherence-free subspace of two qubits against non- 
dissipative errors arising from coupling with unutilized 
(environmental) degrees of freedom. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have developed a prescription for 
calculating density matrices representing one- and two- 
photon states in the presence of unitary frequency- 
polarization coupling. With these methods, wc have in- 
vestigated decoherence in quantum systems arising from 
coupling to unobserved (environmental) degrees of free- 
dom. By inducing frequency-dependent phase errors in 
photon polarization states and tracing over frequency de- 
grees of freedom, we observe information loss as a par- 
ticular form of decoherence. 

In the one photon case, we applied the results of j|] 
and established the feasibility and utility of "bang-bang" 
quantum control of decoherence in an optical system. By 
periodically exchanging the eigenstates of the "environ- 
mental" coupling, photon coherence (polarization infor- 
mation) is preserved. We observed this effect in the case 
in which the strength of the coupling is varied (the "slow- 
flipping" case); the case in which the basis of the cou- 
pling is varied; and the case in which both the basis and 
the strength of the coupling are varied simultaneously. 
Additionally, we demonstrated the utility of bang-bang 
control in the presence of a non-unitary (dissipative) cou- 
pling operation. In all cases, we see a qualitative reduc- 
tion in decoherence for the case in which the quantum 
control procedure is implemented. 

In the two-photon case, we applied general results con- 
cerning the existence of decoherence-free subspaces in 
H to the case of polarization entangled photon pairs 
produced in a down-conversion scheme. Wc found that 
energy conservation in the down-conversion process gives 
rise to frequency entanglement which alters the assump- 
tions underlying the prediction of a DFS. In particular, 
frequency anti-correlation breaks the permutation sym- 
metry of the environmental coupling between opposite 
photon paths by anti-correlating their respective phase 
errors. We found that modifying the experimental appa- 
ratus to restore the symmetry of the phase errors allows 
us to produce a truly collective decoherence process, and 
consequently we recovered the singlet state as a DFS. 

In conclusion, we have investigated two procedures for 
reducing the decoherence of an optical qubit subjected 
to non-dissipative phase errors and non-unitary dissipa- 
tive errors. The first procedure, "bang-bang" quantum 
control, preserves polarization information in a single- 
photon by rapid, periodic control operations which ef- 
fectively average out decoherence effects. The second 
procedure requires entanglement between two photons 
so that the resulting state exhibits symmetry properties 
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which are immune to collective decoherence. These cases 
together demonstrate the feasibility of both active and 
passive protection of coherence properties in an optical 
qubit. 



APPENDIX A: MATHEMATICAL ARGUMENTS 

1. Eigenvalues of U(i) 

Unitarity restricts the possible eigenvalues, Uj(u>,x), 
ofU(x): 

1 = (XilXj) 

= ( Xj \tf(x)XJ(x)\ Xj ) 

= (Xj\Uj{u,x)Uj(u,x)\xj) 

= lUji^x)] 2 . (Al) 

Since Uj(u>,x) is complex with modulus 1, we may write 



(A2) 



where (fj is a real- valued function of lu and x, and the 
subscript refers to polarization mode j. 



2. Operator Relations with R and U 



a. Proof of Equations 11c, 



111 



Sin ce h i is a line ar operator, we need only show that 
Eqs. 11a and lib hold for the basis states, \i. ui) an d 
|X2, oj). From the definitions of U and R (Eqs. pa - 10b ) . 
we have 

BUBlA\xiM = BUB\xxM 
= TW\ X 2,oj) 



(A3a) 



Similarly, 



RWRW|x2,^) = e^ (w) RWR| X 2,( 

= ±e^\ X 2,Lo). (A3b) 

Thus we have, (RU) 2 = ±e zv ^ where the upper (lower) 
sign occurs when R represents a polarization reflection 
(rotation). 
It follows immediately that 



(A4) 



Since equations |lla| and |llq hold for the basis states they 
hold for all linear combinations of basis states, which is 
to say, for all states. 



Derivation of Equations 21a- 22l\ 



All four relations follow e asily from the definiti ons of 
R and Uj (Eqs. |9j]-|t] and paHl8b|) . Equation |la| is 
trivial, since by definition, both operators U\ and U 2 act 
as the identity on |xi). For the other three, we apply the 
definitions in succession. 



U 2 U 1 \ X i) 
U 2 Ui\x2) 

RW2RW1IX1) 



RW2RW1IX2) 



e i(vi(w)+V2(w))|^ 2 ^ 

RW 2 R| X i) 
R^ 2 |x 2 ) 

Re^^| X2 ) 
±e^ (|xi)) 

Ixi) 

RU 2 Re llfil ^\ X2 ) 
±KU 2 e lVl( ^\ X i) 
±Re llpl{uj) \ X i) 



(A5) 
(A6) 



(A7) 



e t(Vi(w)-v>2(w))|^ 2 \ 



(A8) 

Again, if R represents a reflection (rotation) we use the 
upper (lower) sign. From these results, we see that U 2 Ui 
can be treated as a single operator U with an associated 
phase error of <pi(u>) + tp 2 (uj)- Factoring out the global 
phase shift of (p 2 (c j) in the last two relations (indicated by 
the arrow in eqs. A7 and A8), we see that RU 2 RUi can 
be treated as a single operator U QC with an associated 
phase error of (pi(u>) — tp 2 {oS). 



3. Experimental Measures of Coherence 

a. Proof of Equation [j^ 

We begin by defining A in Eq. [2^ to be an equal su- 
perposition of the basis states \xj)' 

1 

72 
1 

Then we have 



|A) = 
|A) = 



(Ixi) 
(Ixi) 



1X2)) 
|X2)). 



(A9a) 
(A9b) 



V{x) = {\\p(x)\\) - {\\ P {x)\\) 

= \{{{XX\+{X2\)P{X){\XX) + \X2)) 

- ((xi|-<x 2 |)^)(lxi)-|x2))} 
= {{xMx)\x2) + (x2Mz)lxi» 
= ((xi|pWIx 2 ) + (xi|pWIx2)*) 

= 2iZe«xi|p(s)|x2» (A10) 
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where Re indicates the real part. 



J 



4. Explicit Calculation of Two-Photon Density Matrices 



We begin by modeling two-photon density matrices including frequency anti-correlation effects. Substituting the 
appropriate phase functions (Eq. |4l"| ) into Eq. we have an explicit expression for p = p(C,C): 



E c v c *ki\x t j)(Xki\ 

i,j,k,l 



de 



e)| 2 K 



e)|VM ; 



'-+e,jC.)+tpj(^-—e,C) — lfk(^ +£,C)-<Pl{~ 



X)) 



}. (All) 



We will now calculate p assuming the initial state matrix elements, CyC^, are given by the density matrix represen- 
tation of the four Bell States, 



|^) = ± \W)), |^) = ± 



(A12) 



In this way, we see the effect of U on each of these states. For brevity, let /(e) = \A(%jf- + e)\ 2 \A(^ — e)| 2 . Then we 
have, 



( i o o ±r> 
oooo 



V±i o o i 



/0 0\ 

110 1 ±10 

2 ±1 1 

\0 0/ 



/ 1 



Vie 1 ^^ 



o o ± e - 4w ° 




1 







±/de/(e)e 

\o 



) 



±/de/(e)e 2 
1 







0/ 



(A13) 



(A14) 



Note that the dependence on e (the variable of integration) drops out of the (11,22) and (22, 11) terms of |<^ ± )(0 ± | 
so that there is no decoherence due to the frequency distribution of the down-converted photons. However, e does 
not drop out of the (12, 21) and (21, 12) terms of |^ ± )('0 ± |. Due to the finite width of the frequency ter m /( e), for 



large £, the integral in these terms approaches zero in the same manner as the Fourier integrals of Sect. II A. Note 
that the results developed above apply only to decoherence in the H/V basis. The states \<f> ) look different in other 
bases, so they are subject to decoherence in general. 

In order to model decoherence in other bases, and also to model non-collective decoherence, we must include rotation 
operators in Eq. All. We define an operator 1Z(0l, 9r) which represents a rotation by 6l in arm L and a rotation by 
R in path R. It can be shown that in the present matrix notation, 



n(o L ,o R ) 



/cos(6l)cos(6r) — sin(^i) cos(0r) — cos(#l) sm(0n) sui(#l) sm(6 R ) 

sin(0[,) cos(0r) cos(0l)cos(0r) —sin" 1 "' '' 

cos(0 L )sm(0 R ) -sm(0 L )sm(O R ) cos(fr L j cos^; 

sm(O L )sin(0 R ) cos(9 L ) sin(9 R ) s'm(0 L ) cos(0 R ) 



\ 



{0 R ) -cos(O L )sm(0 R ) 

I L )COS(0 R ^ ~K\r,(9r)mx(9„) 



-sin(0 L )cos(0fl) 
cos(0 L ) cos(9 R ) 



(A15) 



It follows that 1Z~ 1 (9l,9 r ) = H)(9l,9 r ) = 1Z(—9l,—9 r ). We can now generalize Eq. [59] by writing the reduced 
density operator representing polarization degrees of freedom after an input state p w (0) has passed through crystals 
of thicknesses xl and x R at angles 8l and 9 R in the photon paths labeled L and R: 

p( Xl , X r,6 l ,9r) = 1 7 '6w 3 dw4(Lj 3 \(LJi\n(6 L ,e R )v(x L ,x R )n\e L ,e R ) P ^ 

(A16) 



Eq. A16 can be solved numerically given cpj(u),x) (j G {1,2}), the eigenvalues which define XJ(xl, xr). 

With the same techniques, we can model density matrices including both frequency and path anti-correlation effects. 
Substituting the modified phase functions ( f45| ) into Eq. ^6|, we have an expression which includes both frequency and 
path anti-correlation effects: 



P= E ^ c UXij){Xkl\ j de {|A(^ + e)| 2 |A(^-e)|V( 



X)) 



i,j 7 k,l 



}. (A17) 
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Again, using the four Bell States as input states, we see the effect of U on each of these states. 



( 1 ±1\ 



V±i o o i 

/0 0\ 

1 ±1 

±1 1 

\0 Oy 



1 












2 














{±fdef(e)e- 











/0 







°\ 


1 


1 









2 


ie ! "°^ 


1 









[p 







0/ 



± f de/(e)e 2iei 


1 



(A18) 



(A19) 



Note that in this case, the dependence on e (the variable 
of integration) drops out of the (21, 12) and (12, 21) terms 
of )(V> |- On the other hand, e does not drop out of 
the (11,22) and (22,11) terms of |0 ± )(0 ± |. Again, the 
results developed above apply only to decoherence in the 
H/V basis. The analytical form of p is helpful in devel- 
oping a physical picture of the decoherence process. In 
order to calculate the corresponding result for decoher- 
ence in an a rbitrary basis, we simply set Ot, = Or — 5 
in Eq. A16 and solve numerically. This corresponds 
to the path anti-correlation procedure described earlier. 
These results will be presented elsewhere JO!, however, 
the theoretical predictions are unchanged for the state 



\ip ), since it looks the same in every basis. 
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